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We study the pair production of scalar top-quark partners decaying to a top-quark pair plus large 
missing energy at the LHC, a signature which appears in numerous models that address outstanding 
problems at the TeV-scale. The severe experimental search cuts require a description which combines 
higher-order corrections to both production and decay dynamics for a realistic final state. We do 
this at next-to-leading order in QCD for the first time. We find large, kinematic-dependent QCD 
corrections that differ dramatically depending upon the observable under consideration, potentially 
impacting the search for and interpretation of these states. 



I. INTRODUCTION 

The ATLAS and CMS experiments have recently an- 
nounced the discovery of a Higgs-like boson at the 
LHC H|, thus beginning the experimental exploration 
of the underlying mechanism of electroweak symmetry 
breaking. There are numerous reasons to expect some- 
thing more intricate than the single boson predicted by 
the Standard Model (SM). Naturalness arguments pre- 
dict the existence of new heavy partners of the top quark 
that cancel the quadratically divergent contribution of 
the SM top quark to the Higgs mass, thereby stabilizing 
the hierarchy between the electroweak and Planck scales. 
Additionally, the relic abundance of the dark matter in 
the universe is naturally explained by a stable, neutral 
particle with a mass near the electroweak scale. Many ex- 
tensions of the SM attempt to simultaneously solve both 
of these issues, and contain both a heavy new particle 
with the gauge quantum numbers of the SM top quark, 
and a new discrete symmetry which makes the lightest 
parity-odd particle a good dark matter candidate. Ex- 
amples of such models are the Minimal Supersymmetric 
Standard Model (MSSM) [|| and the Littlest Higgs with 
T-parity [4|. The MSSM contains a spin-0 stop and a 
spin-1/2 neutralino that fulfill the aforementioned roles, 
while the Littlest Higgs with T-parity contains a new 
spin-1/2 fermion denoted as T and a parity-odd partner 
of the photon called an Aq, which respectively serve as 
the top-quark partner and the dark matter candidate. 

A generic production mode of such theories is the 
QCD-initiated pair production of two top partners, fol- 
lowed subsequently by their decay into the SM top quark 
plus the dark matter candidate. This leads to the signa- 
ture 
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where T generically denotes the top-quark partner. Such 
a process could be the dominant signature for supersym- 
metry in 'natural SUSY' models that contain a light stop 
quark and a somewhat heavy gluino Q. Top-quark plus 
missing energy signatures have been considered numer- 
ous times in the theoretical literature 0, 0], and have 
been searched for experimentally [jl-tToj . The current ex- 
perimental limits exclude spin-^ top-quark partners with 
a mass between 300-480 GeV for m Ao = 100 GeV [T(J, 
and scalar partners with a mass in the range 300-450 
GeV for a similar rriA value Q. The proposed theoret- 
ical search strategies, and those utilized experimentally, 
all require an excess in the tail of an energy-related dis- 
tribution, such as ET,miss, the transverse mass of the 
lepton and missing Et if the top-quark pair decays semi- 
leptonically (denoted by My in this manuscript), or the 
effective transverse mass Mr, e // = Ex.miss + Si^r.i, 
where i runs over all observable particles and Et,i = 
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i + Pi t ■ I* been emphasized that variables such 
as Mx.eff may also help distinguish the spin and other 
properties of the top partner [J. 

In this manuscript we wish to improve upon the de- 
scription of the tt + ET,miss signal process, to assist in 
both the search for and eventual interpretation of the un- 
derlying model assuming discovery In all experimental 
searches so far performed, the signal process was modeled 
with leading-order kinematics, and was normalized to an 
inclusive higher-order prediction for stable stops [TJ (3 ■ 
The severe experimental search cuts require a description 
which combines production and decay dynamics for a re- 
alistic final state. An exact next-to-leading order (NLO) 
QCD analysis was not performed, nor were NLO QCD 
corrections considered in the decay of the T particle, or 
in the decay of the top quark. Indeed, such a complex 
new physics signature as that considered here has never 
before been studied with the exact NLO QCD correc- 
tions included consistently through the entire production 
and decay chain. We perform such an analysis for the 
first time for the ti+ Ex^miss signature together with the 
semi-leptonic decay of the top-quark pair. We find large, 
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kinematic-dependent QCD corrections that differ dras- 
tically depending on the observable studied. We show 
that the higher-order corrections behave very differently 
for the Ex^miss, Mt and Mx, e // variables, and identify 
the reasons for these differences. Such effects must be 
accounted for in interpreting the implications of these 
searches. In the following sections we discuss the details 
of our calculational framework, and present illustrative 
numerical results for an 8 TeV LHC. We focus here on 
the scalar top-partner case, leaving a more detailed study 
of both the scalar and a fcrmionic partner to future work. 



II. CALCULATIONAL FRAMEWORK 

We begin by outlining the techniques used to obtain 
our results. We calculate the NLO QCD corrections to 
the process pp — > TT — > bblvjjAaAo by extending the 
framework of Ref. fl3j for top-quark pair production. 
We consider the production of a scalar TT pair which 
is followed by consecutive on-shell decays of T — > tAo, 
t — > bW and W — > Iv/jj. For simplicity, we assume 
that the scalar top partner decays 100% of the time 
through the process T — > tAo. We neglect contribu- 
tions that are paramctrically suppressed by 0(Tt/ttit), 
0(T t /m t ) and 0(Tw /fnw), in each of the decay stages 
respectively. This sequential framework is then system- 
atically promoted to NLO accuracy by calculating QCD 
corrections to the production and decay processes, in- 
cluding all spin correlations in the narrow- width approx- 
imation. If desired, we can systematically improve our 
approximation by allowing off-shell top quarks in the 
decay T — » tA Q — > WbA . We numerically calculate 
virtual corrections for the production process using D- 
dimensional generalized unitarity methods |14j , which we 
extend by adding new tree-level recursion currents in- 
volving scalars, quarks and gluons. Real corrections to 
TT pair production do not exhibit final-state collinear 
singularities and soft singularities are spin- independent, 
allowing us to reuse previous results for top quarks [l3| . 
QCD corrections to the decay T — > A^t are derived an- 
alytically using a traditional Fcynman-diagrammatic ap- 
proach. Similarly to the production process, we can make 
use of existing results for top quarks to treat singulari- 
ties in the real-emission decay process. We subtract the 
soft singularity in T — > A^tg with the dipoles of Ref. [l5[ 
which were developed for the decay t — > Wb retaining a 
finite 6-quark mass. QCD corrections to the remaining 
stages in the decay chain, t —> bW and W — y jj, are taken 
from previous results for top-pair production. To ensure 
the correctness of our calculation we performed several 
cross-checks. First, we confirmed numerically that 1/e- 
poles in dimensional regularization, where e = (4 — D)/2, 
cancel between virtual and real corrections in the produc- 
tion as well as in the decay matrix elements. To check 
the finite parts, the virtual corrections to the process 
qq — > TT have been confirmed by an independent Feyn- 
man diagrammatic calculation for stable scalars T and 



T. Similarly, the virtual correction to the decay process 
T — > A t has been cross checked by a second indepen- 
dent calculation. The implementation of all real correc- 
tions has been checked for independence on the cut-off 
parameter a that controls the resolved phase space of 
the dipole subtraction terms. We also verified the cor- 
rectness of our calculation by numerically comparing to 
the results of Ref. [lj| as implemented in Prospino 2 . 1 
for stable stops in the heavy-gluino limit. We find 
very good agreement for the total hadronic cross section 
at NLO QCD. To further check the implementation of the 
decay stages, we tested factorization properties between 
production and decay matrix elements. This is achieved 
by removing all acceptance cuts on final state-particles 
and integrating over the full phase space. The result is 
compared to a separate evaluation of the product of to- 
tal cross section for stable squarks times their branching 
fraction. We find that the required identities are fulfilled 
within the numerical precision. 

Although we are considering here only a simplified 
model with a single scalar top-partner and a stable spin- 
1/2 particle, we comment briefly on how we expect these 
results to extend to the MSSM. Using Prospino we see 
that this production channel receives negligible gluino 
contribution once its mass exceeds one TeV. Hence our 
results for the production cross section will hold in the 
appropriate parameter region of the MSSM. We have as- 
sumed that the scalar top partner decays entirely through 
T — > tAo, which would be the case in the MSSM for heavy 
charginos. The MSSM contains two stop states, the part- 
ners of the left-handed and right-handed top quarks, that 
mix to form the physical eigenstates. The mixing affects 
the couplings relevant for the T — > tAo decay, which are 
free parameters anyway in our study. The contribution 
from the second stop eigenstate would need to be added 
in our study only if its mass is near the lightest eigen- 
state. We therefore expect our conclusions to hold for 
stop production in the MSSM also. 



III. NUMERICAL RESULTS 

To illustrate the impact of the higher-order QCD cor- 
rections, we present results here for several distributions 
at an 8 TeV LHC for pp -> TTA A ->■ bblvjjA A . 
For each observable we show three predictions: the LO 
result, the full NLO result with the QCD corrections im- 
plemented throughout the entire decay chain, and the 
result with NLO corrections included in the pp —> TT 
production process only. Motivated by the experimental 
cuts in the ATLAS analysis @, , we have applied the 
following cuts: 

ARj = 0.4, PTb > 30GeV, 

\y b \ < 2.5, p Tj > 30GeV, 

|y 3 - 1 < 2.5, pt/>20GcV, M<2.5, 

E T>mi ss > 150 GeV, M T > 120 GeV, (2) 
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where Mt is defined as 

M T = 2pTiE T , miS s (1 - cos(A0)), 

with A(j> being the azimuthal angle between the lep- 
ton and missing energy vector. We have used the 
MSTW2008 parton distribution set [3 with the cor- 
responding choices of the strong coupling constant: 
a SLO (M z ) = 0.13939 and a SNLO (M z ) = 0.12018, which 
are subsequently evolved to the scale choice \x using 1- 
loop and 2-loop running at LO and NLO respectively. 
We present results for the following choices of the mass of 
the top partner and the stable particle A : (mr, rnA ) = 
(500 GeV, 100 GcV). The scale fi for the central value 
in each distribution is taken to be 500GeV. We note 
that we have also studied the mass point (mr,W4 ) = 
(250 GeV, 50 GeV) which belongs to the compressed spec- 
trum region that is currently not excluded, and have ob- 
tained similar results to those presented here. Our choice 
of the left- and right-handed couplings of the top partner 
to the top and A$ are as follows: 

gn = cr m t /v; g L = c L m t /v, (3) 

with c R = 3/10, c L = 1/10, m t = 172 GeV and v = 246 
GeV. We note that this choice does not have a strong 
theoretical motivation and is only meant to illustrate the 
impact of higher-order effects. 

We begin by showing the distributions for the two pri- 
mary variables used in a recent ATLAS search for this sig- 
nature in the semi-leptonic mode the missing trans- 
verse energy distribution in Fig. [TJ and the transverse 
mass of the lepton and missing energy in Fig. [5J Several 
features are apparent from these plots. First, the cor- 
rections are generally large, with the ratio of NLO over 
LO (defined as the if -factor in the plots) ranging from 
approximately 1.6 to 2. Second, neglecting the QCD cor- 
rections in the decay of the TT final state leads to an 
overestimate of the size of higher-order corrections by 
about 10%. Finally, we note that the kinematic depen- 
dence of the corrections differs for the two observables. 
For ET,miss, the if- factor begins at around 1.8, plateaus 
near 2 for Ex^miss ~ 600 — 800 GcV, and gradually re- 
duces for higher values. For the transverse mass, the 
if- factor is 2 for Mt near the lower kinematic limit. It 
then monotonically decreases to 1.5 for Mt ~ 800 GeV. 

The size of higher-order corrections and kinematic de- 
pendence can differ even more dramatically for other ob- 
servables. We demonstrate this in Fig. [3] by showing the 
effective transverse mass of the entire final state. The 
shape of this distribution for large Mx.eff has been sug- 
gested as a useful diagnostic tool to determine the top- 
partner spin Q. The if- factor rises from below unity to 
values of 3 or more at large M^e//. We can qualita- 
tively understand this behavior as follows. Our analysis, 
and the one performed by ATLAS, demands four jets in 
the final state. At high M^e//, the top quarks become 
boosted, and the probability of finding four jets at LO 
with an anti-fey separation parameter of Ai? = 0.4 is re- 
duced. This issue is alleviated at NLO by the presence 
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FIG. 1: The missing energy at LO and at NLO for 
(mT,mA a ) — (500GeV, lOOGeV). The upper panel shows 
the distributions while the lower panel shows the if -factors, 
defined as the ratio of NLO over LO, both with and without 
corrections throughout the entire decay chain. 
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FIG. 2: The transverse mass of the lepton and missing energy 
at LO and at NLO for (m T , m Ao ) = (500 GeV, 100 GeV). The 
upper panel shows the distributions while the lower panel 
shows the if -factors, defined as the ratio of NLO over LO, 
both with and without corrections throughout the entire de- 
cay chain. 

of additional radiation. This explanation is confirmed by 
checking that the if -factor decreases when the anti-fey 
parameter is reduced to Ai? = 0.1. A LO calculation of 
the decay process fails to properly describe this observ- 
able. It may be possible to avoid such large QCD cor- 
rections by incorporating techniques for tagging boosted 
objects into the search strategy. We also study the inte- 
grated cross sections and acceptances obtained by using 
the cuts in Eq. ([2]) at LO and NLO. Our results are shown 
in Table |U The acceptance increases from 0.20 to 0.29, 
i.e. by 45%, when going from LO to NLO. This would 
shift the excluded cross section determined by the exper- 
imental collaborations by a similar amount, significantly 
impacting the interpretation of experimental results. We 
note that the scale variation of the acceptance almost 
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Mr.eff [GeV] 

FIG. 3: The effective transverse mass of the final-state sys- 
tem at LO and at NLO for (rn T ,m Ao ) = (500 GeV, 100 GeV). 
The upper panel shows the distributions while the lower panel 
shows the if-factors, defined as the ratio of NLO over LO, 
both with and without corrections throughout the entire de- 
cay chain. 



Cross section 


no cuts 


with cuts 


acceptance 


o"lo 


4.571^ fb 


0.9lt° ;f fb 


0.201° 


CNLO 


6.071°;®® fb 


1.77«;f 7 fb 


29+ 02 



TABLE I: The cross section and acceptance values for 
(rn T ,mA ) = (500 GeV, 100 GeV) at LO and NLO using 
the cuts defined in Eq. ^ for ^fs = 8TeV. The central 
value corresponds to fiR = /j,f = OTt and the upper and 
lower uncertainties correspond to = fj,F = 2m,T and 
= A*f = mr/2 respectively. 



completely vanishes at LO, and does not accurately re- 
flect the impact of higher-order effects. The scale depen- 
dence of the inclusive cross section decreases from 32% at 
LO to 14% at NLO. However, once the cuts are included 
the reduction is only from 33% to 23%. Determining the 
uncertainty on the cross section using stable scalars leads 
to an underestimate of the theoretical error. 



IV. CONCLUSIONS 

In this manuscript we have analyzed the production 
of scalar top-partners and their subsequent decay to a 
top-quark pair plus large missing energy at the LHC, a 
signature which appears in numerous models that ad- 
dress outstanding problems at the TeV-scale. We have 
included the NLO QCD corrections throughout the en- 
tire production and decay chains. The impact of higher- 
order corrections depends strongly on the observable un- 
der consideration; the differential if -factor for various 
observables can differ by more than a factor of two in rele- 
vant phase space regions. Current experimental analyses 
which include higher-order corrections only by normaliz- 
ing to inclusive results could be significantly misidentify- 
ing which regions of parameter space are exlcudcd. We 
encourage the experimental collaborations to reconsider 
the allowed parameter space regions in light of these re- 
sults. In future work we plan to thoroughly study both 
scalar and fermionic top-quark partners, and determine 
how NLO QCD affects the discrimination between these 
two possibilities. 
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